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Abstract 

While it is well known that complete positivity guarantees the fulfilment of the 
second law of thermodynamics, its possible violations have never been proposed as 
a check of the complete positivity of a given open quantum dynamics. We hereby 
consider an open quantum micro-circuit, effectively describable as a two-level open 
quantum system, whose asymptotic current might be experimentally accessible. This 
latter could indeed be used to discriminate between its possible non-completely positive 
Redfield dynamics and a completely positive one obtained by standard weak-coupling 
limit techniques, at the same time verifying the fate of the second law of thermodynamics 
in such a context. 


1 Introduction 

The so-called open quantum system paradigm is a remarkably successful way to cope 
with quantum systems weakly interacting with their environments; it has been applied to 
atomic and molecular physics, quantum optics, quantum chemistry and condensed matter 
physics [U El |3l S] . 

In most physical applications there are no initial statistical correlations between system 
and environment and their interactions are weak; in such cases, a reduced dynamics of 
semigroup type for the open quantum system alone can be derived by tracing away the 
environment degrees of freedom, through the so-called weak-coupling limit techniques [T]. 
The corresponding time-evolution is irreversible and characterized by dissipation and noise; 
furthermore, it does not only preserve the positivity of the time-evolving states of the 
system, but is also completely positive laEiii!- Positivity preservation ensures that the 
eigenvalues of the density matrices describing the states of the open system remain positive 
in the course of time so that they can be interpreted as probabilities. Instead, complete 
positivity guarantees the positivity at all times of any entangled state of compound systems 
consisting of the open quantum system of interest and any dynamically inert finite level 
system to which the former may happen to be statistically coupled. 
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Notice that, when initially there are no correlations between system and environment, 
the reduced dynamics that results from eliminating the environment degrees of freedom au¬ 
tomatically consists of completely positive maps; only, the latter do not in general compose 
as a semigroup. In order to obtain a forward-in-time composition law, a desirable physi¬ 
cal property in a weak coupling regime, one usually performs various kinds of Markovian 
approximations that often lead to loss of complete positivity and even of positivity [8]. 

The justification for asking that the semigroup dynamics be completely positive because 
of a general and uncontrollable correlations between the open system and arbitrary inert 
ancillas, is often criticised in the literature as an abstract mathematical artifact u uni 
mm- Indeed, because of complete positivity, the generator of the reduced dynamics 
assumes the so-called Lindblad form to which there correspond specific physical constraints: 
typically, a hierarchy among the decay times of the various entries of the density matrices 
that describe the state of the open quantum system Without complete positivity, 

less constrained dynamics emerge that render easier the occurrence of so-called “quantum 
miracles”, like, for instance, the beating of classical bounds in the efficiencies of energy 
transport processes. 

So far, the completely positive character of a dissipative dynamics has only been ver¬ 
ified by checking, in very few cases, the above mentioned hierarchy, in general a difficult 
experimental task. Instead, in the following, we offer a different strategy based upon the 
thermodynamic behaviour of open quantum systems. Indeed, complete positivity was soon 
recognised to imply the positivity of the internal entropy production as required by the 
second law of thermodynamics naiisiin!; although there is no proof that complete pos¬ 
itivity is also necessary to avoid conflicts with thermodynamics, nevertheless one expects 
that, in absence of complete positivity, the second law of thermodynamics might indeed be 
violated. However, no instances either theoretical or experimental have so far been inves¬ 
tigated in this respect. In the following, we propose a concrete physical context where to 
study a possible conflict between a non-completely positive dynamics and the second law 
of thermodynamics Q 

We shall focus upon a model consisting of three electrons moving in a micro-circuit 
consisting of a three site loop, under the action of a periodical driving and of a weak coupling 
to a thermal bath B modelled as a collection of free harmonic oscillators in equilibrium at 
temperature T. This system, effectively describable as an open 2-level system, has been 
introduced to study the dependence on the driving of the asymptotic current that sets in 
because of the thermal environment m- In the following, we shall study its behaviour 
from a thermodynamic point of view and show that the non-completely positive Redfield 
dynamics considered in m to model the behaviour of the system, violates the second law 
of thermodynamics. On the other hand, by adapting standard weak-coupling techniques [3 
mm. originally developed for time-independent system-environment Hamiltonians, to the 
present driven system, we derive a completely positive reduced dynamics with explicitly 
time-dependent generator which gives rise to temporal patterns of the current supported 
by the microcircuit that significantly differ from those presented in m- The violations 
of the second law of thermodynamics by the Redfield-type time-evolutions can then be 
measured by experimentally accessing the current time-behaviour, thus offering both a test 

preliminary investigation on these topics has been reported in m- 
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of complete positivity and of the fulfilment of the second law of thermodynamics in open 
quantum system dynamics. 


2 Thermodynamics of open quantum systems 

In this section we provide a brief overview of the approach to the thermodynamics of open 
quantum systems as formulated in m- 

A quantum thermodynamical system is taken to be a finite n-level quantum system S 
subjected to a periodically driven time-dependent Hamiltonian H{t) accounting for work 
being performed on 5 in a cycle. On the other hand, heat flows into and out of S because 
of weak interactions with an environment B whose effects are supposed to be effectively 
described by a master equation with explicitly time-dependent generator Lt, 

^ = L,|e(t)l = -i Wit), (.{()] + D,[£.(*)] . (1) 

where Q{t) is the density matrix describing the state of S and h has been set equal to 1. 

The dissipation and noise contributed by the heat bath B are accounted for by the 
time-dependent Lindblad-like contribution 

j,k 

with suitable nxn matrices Vj. The master equation gives rise to a two-parameter semigroup 
of trace-preserving maps 

7m t Tr(7i,^£») = Tr£» , t>s>0, 

where T denotes time-ordering. The complete positivity of such maps is guaranteed by the 
Kossakowski matrix, Kt = [A'jfc(t)], being positive-definite for all t > 0 [1]. 

Equipped with the formalism of open quantum systems, one can formulate the quantum 
versions of the laws of classical thermodynamics as follows. 

• Zero-th law of thermodynamics: it regards the fact that systems initially at different 
temperatures eventually reach thermal equilibrium. In the usual open quantum system 
setting, namely when there is no explicit time-dependence in the generator of ([l]); Ibj — L 
and H{t) = H, the zero-th law corresponds to the system S reaching the thermal equilibrium 
state 

e-/3J^ 1 

= Tr(e-/3^^) ’ ^ ^ ^ ^ ° ’ 

at the temperature T of its heat bath environment, k being the Boltzmann constant. In the 
present thermodynamic setting with an explicitly time-dependent generator Lt, the zero¬ 
th law amounts to the request that the instantaneous Gibbs state of S at time t (at the 
reservoir temperature), Qp{t), be a steady state with respect to the reduced dynamics: 

e-pHit) 

= ’ Li(t)] = 0 . (3) 
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This request is physically justified when the system-bath coupling is weak and the external 
driving by H{t) very slow in comparison with the dissipative time-scale. 


• First law of thermodynamics: it concerns the relation between the rate of work done 
by the system, * , and the rate of heat absorbed from the bath, with the internal 


dt 

energy variation (energy balance): 


dt 


dEt _ dQt dWt 
dt dt dt 


where Et = Tr {^g{t)H(t)). The work per unit time done by the open quantum system due 
to the time-dependence of the system Hamiltonian is 


dWt 

dt 


-Tr 



Then, the heat exchanged per unit time by the system amounts to 


( 4 ) 


^ ^ d{Et + Wt) ^ Trr mti . 

dt dt \ dt y 

By means of the master equation 0 . this expression can be recast as 
^ = Tr(F(t)Li[^,(t)]) = . 


( 5 ) 

( 6 ) 


• Second law of thermodynamics: it regards the variation of the internal entropy of a 
system, namely the entropy variation which is not due to heat exchanges between the 

SQ 

system and the environment and thus not of the form The second law asserts that the 
internal entropy cannot decrease in time. 

In order to cast this request in a mathematical form, firstly one introduces the total 
entropy of the system as given by the von Neumann entropy of its state Q{t), 

S{Q{t)) =-kTi(^Q{ t) log g{t)j . (7) 


Then, one considers that, as outlined above, its variation in time, 


dSjgjt)) _ d5int(t) d5ext(^) 

dt dt dt 


( 8 ) 


is due to two terms: an external entropy rate related to the heat exchange with the envi¬ 
ronment, 


d5ext(t) 1 dQt 


dt 


T dt 


Kl3TT(^H{t)ht[g{t)]^ , 


(9) 


and an internal entropy rate 


dSint(t) 

dt 


that we shall denote by a{g{t)) for sake of simplicity. 
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Definition 1. The internal entropy production in a quantum thermodynamic system evolv¬ 
ing in time according to the master equation Q is given by 

cr{e{t)) = ~ = - kTt (ht[Q{t)] (^log Q{t) + (3 . (10) 

The second law of thermodynamics amounts to the request that cr(^(t)) > 0 for all t > 0. 

The positivity of the internal entropy production is guaranteed by the complete posi¬ 
tivity of the maps generated by the master equation 0. This can be seen as follows: using 
the definition Q of the instantaneous Gibbs state, one has 

TT(ht[g{t)]PH{t)^ =-Tr(^L4[^5(t)] log^^(t)) ; 

then, one can rewrite 

o'(^>(t)) =-KTr^L4[^(t)](^log^>(t) - log^;3(t)^^ . (11) 

The sign of a{g{t)) is accessed through the relative entropy of two density matrices ^> 1 ^ 2 , 

*S’(£'i|£' 2) := Tr^£>i(^log£»i - log£»2)^ , (12) 

for it monotonically decreases under completely positive, trace-preserving maps |20j . 

Indeed, for each fixed t > 0, one can use the generator L* of the master equation 0 
and construct a semigroup of maps = exp (sL^), s > 0. Due to the Linblad form Q of 
Lt the maps are completely positive and trace-preserving; moreover, from L,t[g^{t)] = 0 
it follows that As[^_( 3 (t)] = gjsit) for all s > 0. Then, 

S{As+s[g{t)]\g^it)) = 5(As+5[£)(t)]|As+5[^;3(t)]) = S{As o A^[£i(t)]|A5 o A5[£)/3(t)]) 

< S(A4£.(t)]|A4£.^(t)]) = 5(A^[^(t)]|£»/3(t)) . 

By taking the derivative of —S'(As[^j(t)]|^)^(t)) with respect to s at s = 0, it follows that 
0 < -/t d5(A^[g^tj]|g^(t)]) =-KTT(ht[git)](^logg{t)-\oggi 3 {t)^'^ = a{g{t)) . (13) 

For the interpretation of the term —Tr^)(t) log ^//^(t) in the relative entropy as a heat 
flow, it is crucial that gjsit) be a Gibbs state 0 such that ht[gj 3 {t)] = 0. In many cases of 
physical interest, however, the steady states gstit), Lt[^Jst(t)] = 0, need not be Gibbs states. 
In these cases, the notion of entropy production can be extended following the argument 

of US]. 

Definition 2. In the case of a quantum thermodynamical system undergoing a dissipative 
dynamics with a stationary state gst{t) which is not a Gibbs thermal state, the expression 
of the internal entropy production in Definition is generalized to 

a{g{t)) := -k ^S{g{t)\gstit)) = -«;Tr|^Lt[£)(t)](log£>(t) - log£>st(t))^ , (14) 

with aig{t)) > 0 for all t > 0 expressing the second law of thermodynamics. 
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That the expression (14) can be interpreted as a bona fide internal entropy production 
is based upon its vanishing when the dynamics is unitary and its convexity as required by 
thermodynamical stability |16) . 


3 Open quantum micro-circuit 


In this section, we study the second law of thermodynamics in the following solid state 
model introduced in |19j : a triangular loop micro-circuit comprising three equal, single 
orbital quantum dots with three electrons that can jump from one to the other under the 
action of a periodic modulation of the transmission amplitudes and gate voltages. The 
purpose of m was the study of the asymptotic current which sets in when the micro- 
circuit S is placed in weak interaction with a thermal bath B consisting of non-interacting 
harmonic oscillators. We briefly present the main characteristics of the model indicating 
how it can be reduced to the study of an open two-level system. 

Let a, b, c label the three quantum dots; in absence of thermal bath, their dynamics is 
generated by a Hamiltonian that, in the standard basis |o), \b) and |c) is chosen of the form 


Ho{t) 


'£a{t) 

-70 

-70 

-70 

£b{t) 

-70 

,-70 

-70 

£c(i) 


where 70 > 0 is the hopping amplitude from one dot to the other and £a,b,c{t) are the 
following external biases: 

£a{t) = —A cos(Ht) , £b{t) = —A cos (ilt — 2 tt/ 3 ^ , edt) = —A cos + 27r/3^ 

with A > 0 such that A <C 70 and H the external driving frequency. The unperturbed 
Hamiltonian Hq obtained from Hq^) by setting the biases equal to 0 has a ground state 


| 0 ) = 


|a) T |&) “ 1 “ |c) 


with energy —270 and the following two degenerate orthogonal excited 


states with energy 70 


x) 


1^) - |c) 

v/2 


\y) 


2\a) - \b) - |c) 

v/6 


In the hypothesis of neglecting electron-electron correlations, when three electrons move 
over the micro-circuit, two of them are expected to occupy the ground state, while the state 
of the third one belongs to the orthogonal subspace. 

As the biases are considered as small perturbations, one may neglect transitions off 
the ground state and thus restrict to considering the third electron described by the two- 
dimensional subspace spanned by |x) and \y). The restriction of Ho(t) to such a subspace 
yields a pseudo-spin Hamiltonian 

H(t) = — ((T 3 cos Hi-|-fji sin , (15) 

where o'i^ 2,3 are the standard Pauli matrices with respect to the basis |x) and \y). 
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The description of the micro-circuit can thus be effectively reduced to that of a two- 
level system with time-dependent Hamiltonian. Within this setting, the electronic current 
at time t sustained by the micro-circuit is measured by the mean value of (T 2 with respect 
to the time-evolving two-level system density matrix g{t): 

I{t) = IoTtsi0{t)(r2) , /o = ^ , (16) 

where e is the electron charge. Indeed, the current is induced by the hoppings of the 
electrons between adjacent dots and thus it is described by the self-adjoint matrix 

/ = ie'yo(\a){b\ - |6)(a|^ . 

By rewriting I in the orthonormal basis |0), jx), \y), and neglecting transition terms con¬ 
necting the unperturbed ground state |0), and the excited states |x) and \y), one finds 

The effects of non-negligible interactions of the micro-circuit with the environment that 
surrounds it will be described by coupling the effective degrees of freedom to a heat bath B 
consisting of independent harmonic oscillators at temperature T. Namely, one considers a 
total Hamiltonian of the system S + B oi the form = H{t) + Hq -|- Hsb, where 



is the bath Hamiltonian with and ^ = 1)3, position and momentum operators of 
the bath oscillators. These are in turn coupled to the micro-circuit degrees of freedom by a 
spin-Bose interaction Hamiltonian: 

A Hsb = A n\/2mUJn ® q^,n , 

?= 1,3 n 

where A is a dimensionless coupling, while the constants A^ are suitable energies associated 
with the bath spectral density. 

Remark 1. Here it is important to notice that the system-bath coupling is taken to be 
homogeneous, i.e. the coupling constants A^ do not depend on the index ^ of ci^. 

The states qsb of the compound system S + B will evolve in time according to the explicitly 
time-dependent Liouville-von Neumann equation 

=-i [H{t) + Hb + XHsb , QSB{t)] ■ (17) 

Notice that the system S Hamiltonian can be recast in the form 

H{t) = ^R{t)a3 R\t) , R{t) = _ 
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Then, by going to a rotating frame by means of the unitary matrix R{t)^ the time dependence 
in 0 can be moved from H{t) to the interaction term Indeed, by setting 

= R\t) Qssit) R{t) , (18) 


its time-evolution equation reads: 

dgsBit) 


dt 


= -i [HeS+ Hs XHsB{t), 


with a new time-independent system Hamiltonian 

Afj3 - Qa2 
Hes= -^- 

and explicitly time-dependent interaction term 

HsB{t) = EE \/2rnuJ^cr^{t) (g) = R"' {t)a^ R{t) 

?=i,3 n 


(19) 


( 20 ) 


( 21 ) 


Remark 2. In the standard open quantum system approach, with no explicitly time- 


dependent global Hamiltonian, an equation like (19) is the starting point for the applications 


of the so-called weak-coupling limit techniques [T] that lead to master equations that involve 
only the degrees of freedom of the system S . These techniques are based on the assumption 
that the coupling constant is small, A <C 1, and that the initial state of the compound system 
S + B has the factorized form gj^, with g^ a stationary state of the environment. Then, 
on a slow time-scale r = the true dynamics of the open quantum system can be well 
approximated by a semigroup of trace-preserving maps 7 * acting on the density matrices of 
S. Only if performed with due care, the approximated reduced dynamics results completely 
positive; in most cases, too rough manipulations yield maps that do not even preserve the 
positivity of states. In [ 8 ], the source of the problem is identified in the too fast oscillations 
in time due to the system Hamiltonian dynamics. It is also shown there how a suitable 
time-average recovers complete positivity. In Appendix A and B we review these issues 


adapting the arguments of [ 8 ] to the explicitly time-dependent equation (19). 


An instance of a reduced dynamics which fails to be completely positive is provided by 
the one adopted in |19j to model the time-evolution of the three dot micro-circuit. There, 


.-PHb 


is the 


starting with an initial uncorrelated initial state g ® ga, where gn = 

^ Tr(e"'^"«) 

thermal bath equilibrium state at temperature T, a Redfield-type master equation for the 
micro-system alone is presented: 


dg(t) 

dt 





r-l-oo f 

- dn G(u) 

, Ues{u) cr^{-u) g{t) 



- 


+ G*{u) 


g{t) C/eff(tt) a^{-u) Ul^{u ), 


( 22 ) 


In the above expression, g{t) = Trg^^) 5 g(t)j is the density matrix describing the state of 
the micro-circuit S at time t, 


Ues{u) = exp {-i u Reff) 


(23) 















is the unitary time-evolution due to the system S Hamiltonian while the complex 

function 

, , , w/3 . . , , 

cos {Louj coth —- I sm [tdu) 


G{u) = / da;J(w) 
Jo 


is constructed with the bath two-point correlation functions 

cos(a;„jtt) coth — isin(w, 


Trg ( 6 ° Q^2,ri2 6 


and the bath spectral density 


ni ‘ 


2muin 


= Xl5{u}-UJn) . 


(24) 


, (25) 


(26) 


For an infinite heat bath, J{uj) is taken to be of Ohmic form with cut-off frequency lOc- 

J(lo) = uj ex.p{—uj/iOc) ■ (27) 


Remark 3. Notice that, although the Liouville-von Neumann equation (19) depends ex¬ 


plicitly on time, the master equation (22) resulting from the Markovian approximation does 


not. As shown in Appendix A, this is the consequence of the homogeneous coupling between 
system and bath (see Remark]^. 

Because of the disappearance of the explicit time-dependence, the reduced dynamics 
is represented by a one-parameter semigroup {'yt}t>o of trace-preserving maps 'jt which, 
however, are not completely positive. 


Proposition 1. The semigroup generated by (22) consists of non-completely positive maps. 


Proof This is easily checked by setting 

/ + 0O 

dr G(t) Wi (r) , Wi (r) = t 4 ff(r) fii (-r) C 4 g(r) 


r'+OO 


= 


drCir) lF 2 (r) , W 2 {t) = t 4 ff(r) C 4 g(r) 


Then, the right hand side of (22) can be recast in the form 

h^-^[g{t)] = -i \h,s + , ^,(t)l + o^^<^[eit)] 


with a Lamb-shift Hamiltonian term 


Hff = ^ - ^3' ^4 


1 


Ai/t 


Ai/t 


and a purely dissipative contribution 

4 


E Kjvkg{t)vJ - l{vJVk,g{t)} 


j,k=l 

with Kossakowski matrix given by iL = [Kjk] = = 


0 aij’ " VI 0. 
is not positive definite then the generated maps is not completely positive [5]. 


V 2 = ai (28) 
H 4 = as . (29) 

(30) 

(31) 

(32) 

Since K = [Kjk] 
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3.1 Failure of positivity preservation by the Redfield dynamics 


As previously mentioned, reduced dynamics of Redfield type fail in general to preserve even 
the positivity of the time-evolving density matrices: this is indeed the case also for the 
time-evolution used in m- The appearance of negative eigenvalues in the spectrum of Q{t) 
is typically occurring at short times. Indeed, given the master equation (22) with generator 
in the form (30), an expansion of the solution p(t) to first order in t starting with a pure 
state |V’)(t^| yields: 

Q{t) ~ |'0)(V^| -h . 

By choosing \(p) £ such that {(pltp) = 0, one finds that only the first term in (32) 
contributes to the mean value of Q{t) at short times: 




A 


^,0 


= TZe 




(33) 

(34) 


If an orthonormal basis {ip), \(p) can be found in such that < 0, then, for small times, 
the projector \'ip){'ijj\ transforms into an operator of trace one with one negative eigenvalue, 
otherwise all mean values (33) should be non-negative. 


Choosing 


|0) -|- Q:| 1) 

\/l + 


and 


a 


•\/l + 


a\ 


, where a G C and ctsIO) = |0), 


o'sll) = -|1), one gets where 

A(a) = 7^e{((a*)2-l)(a((0|Fl|0)-(l|Fl|l)) + a2^0|Fl|l)-(l|Fl|0)) 


+ 2a* a((0|F3|0) - (1|F3|1)) + a"(0|F3|l) - 


■ LRed 


at short times, A (a) must be non- 


In order not to violate the positivity of e* 
negative for all a. 

Setting C = cos(nr), S = sin(nr), c = cos(a;efrr) and s = sin(a;efrr), where Wefr = 
-|- A^, an explicit evaluation of the 2x2 matrices ITi ,3 in (|28|) and (|29l) yields 


/ 


IFi = 


QlCs-^ 

^eff 


‘^eff 


SA^ 


cC-i 


A Cs- 


=sn' 


cC -|- i 


SAP 


+ 


A Cs- 


csn' 


^eff 


I ssri 

tdeft 


IFs = 


ca 2 o(sS-r^) 
- ^ - 


^eff 


V 


cS -|- i 


■'"(•S+S) CAS! 


cS — i 


^eff 




CsO 

^eff 



‘A’eff 

SA2 

Cs-^) 

M 

dJeff 

A{.s+fs) 

CAn ) 

‘A’eff 


CA2 ^ 

7-7T 

(sS+^) 




ssn 


\ 


CsQ 

‘A’eff 


/ 


By means of these expressions the matrix in (28) and (29) can be numerically computed 
as well as the behaviour of A (a). The following figure shows the complex a’s for which 
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A (a) < 0 and thus the amount of pure states that are sent out of the Bloch sphere at small 
times by getting a negative eigenvalue. 


AW 

(W’+l)^ 



Figure 1: Presence of negative eigenvalues at short times: A(a) < 0. 


3.2 Completely positive dynamics 

A fully physically consistent time-evolution for the reduced system can be obtained by 
means of a more careful manipulation through standard methods in open quantum system 
theory, thus avoiding drawbacks as absence of positivity and of complete positivity. These 
methods, known as weak-coupling limit, are essentially based upon eliminating too fast 
oscillations due to the dynamics generated by by means of a time-average; however, 
they are developed for time-independent system-environment Hamiltonians, whereas in the 
present case we deal with a driven two-level system whence the standard weak-coupling 
limit techniques must be adapted to the new setting. This is done in Appendix B, where it 
is shown in detail how to derive a master equation of the form 


dg(t) 

dt 


= L[^)(t)] = -i + X^Hls , Q{t) + A2D[£)(t)] 


(35) 


that generates completely positive maps. 

In order to compare the previous master equation with the one in (22), it proves conve¬ 
nient to recast them both in a vectorial form using the so-called Bloch-representation with 
respect to a new triple of Pauli matrices 


CJl = Ul 


0-2 = 


Acj 2 -I- fids 

^eff 


0-3 = 


^efF . 


Ads - ^<72 


^eff 


(36) 


so that the effective Hamiltonian (20) becomes i^eff = “^ds. 


One can then expand the time-evolving density matrix as follows, 

6it) = + +^2(i)+ ^’3(^) < 1 : 


(37) 


2=1 
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and identifies it with the real Bloch vector |^?(^)) with components (1,ri(t),r2(t),r3(t)). 

Then, the action of the generator in (30) can be rewritten in terms of a 4 x 4 

matrix pRed that 


(38) 



/ 

1 \ 



/ 1 

\ 

/ 0 

0 

0 



d 

ri(t) 

= -2C 

Red 

ri(t) 

^Red _ 

Tio 

Til 

Ti2 

Ti3 


dt 

r2(t) 


r2{t) 


T21 

T22 

T23 








V-^30 

T31 

T32 

T33/ 


where the contributions 











/o 

0 

0 

o\ 


^0 

0 

0 

0 


U rr - 


0 

0 

Wefr/2 0 

-l/Red 

0 

0 

T-Li2 

Hi3 



0 


0 

0 

) — 

0 - 

'H 12 

0 

^23 



VO 

0 

0 

0/ 


vO - 

T~il3 

—'^23 

0 

/ 


(39) 


come from the two commutators in (30), while the dissipative term corresponds to 


-j^Red _ 


/ 0 0 0 0 \ 

/Cio /Cii /Ci2 /Ci3 

/C20 ^21 A^22 ^23 
V^30 ^31 A^32 ^33/ 


(40) 


The coefficients of these matrices are explicitly given in Appendix A. Instead, in Appendix 


B, the generator of the master equation (35) is proved to correspond to the 4x4 matrix 


C = Heff + X^I-Lls + where Hes is as in (p^, while 


T-Lls = 


/o 

0 

0 



( ° 

0 

0 


0 

0 

T-Li2 

0 

V = 

0 

/Cii + /C22 

0 

0 

0 

—T~il2 

0 

0 

1 ^ 

0 

0 

^11 + ^22 

0 

Vo 

0 

0 

0^ 


V ^30 

0 

0 

^33/ 


(41) 


3.3 Stationary states and cnrrent 

The vectorial form of the master equation in the Bloch representation is particularly suited 
to finding the stationary states Qst such that L[^3st] = 0; indeed, they correspond to the 
Bloch vectors (ro,ri,r2,rs) that solve the linear system 


Cjfiri, = 0 , z = l,2,3, ro = 1 . 

/^=0 


The unique stationary state of the master equation (35) is easily found using (41): it 
corresponds to the Bloch vector 


^st = 


1 + rfas 


r®* — — 
^3 — 


^30 

^33 


(42) 


Therefore, in the case of the completely positive reduced dynamics generated by (35), all 
initial states of the open quantum system will asymptotically tend to ^st- 
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On the contrary, in the case of the non-completely positive Redfield dynamics generated 
by ( 22 ), the stationary state will not in general be unique and the manifold of asymptotic 
states will depend on the explicit form of the coefficients in the 4x4 matrices and 
In |19] . based on numerical evidences, it is argued that, with the following choices of coupling 
constant, A = 0.005, bath temperature, T ~ 0.006 K, pumping amplitude, A ~ 8 GHz, 
pumping frequency, H, and Ohmic cut-off frequency, Uc, such that Wc/A = 10^, H/A = 2 , 
the stationary state in (42) can be taken as an approximation of the true stationary state 
of the Redfield type master equation (22). However, as we shall show below (see Figures]^ 
and 1 ^ , the asymptotic behaviours of the two dynamics give rise to different values for the 
asymptotic current supported by them. 


In the non-rotating frame, the stationary states become time-dependent steady states. 
Indeed, the physical reduced dynamics that sends the initial state g into g{t) is obtained 
from the one sending g into g{t) by going back to the non-rotating representation, where 
g{t) = R{t) g{t) R^t). However, since R{t) leaves a 2 invariant, the currents supported by 
the micro-circuit in the states g{t) and g{t) are the same: 


It = IoTr\^g{t)a2j = IoTr\^g{t) a2j ■ (43) 

In particular, the completely positive time-evolution predicts an asymptotic current sup¬ 


ported by the stationary state which, using the explicit expressions (|87[ ) and (94) of the 
coefficients /C 30 and /C 33 computed with the Ohmic spectral density (|27l), reads 


1st — In 


where, 


/C30 ^30 _ (^eff ~ J-I- -k (cUeff -k J- 

^eff ^33 ^33 (^eff H)^C-|_./-|_ -\- (cUeff -\- H)^C_(/_ 

( /3(Weff ± H) 


(44) 


J± = J{oJett i H), c± = coth 




Asymptotic and thermal states In the expressions ( |11[ ) and (14) for the internal en¬ 
tropy production there appears either a Gibbs state associated to a heat flux and to the 
external entropy production, or a stationary (steady) state. 

In the thermodynamic analysis of the present model, we have three asymptotic states: 


the stationary state gst in (4:2), the Gibbs state associated with the micro-circuit Hamilto¬ 
nian Heff, 


efj ^ e ^ 1 - tanh (a;eff^/ 2 )d 3 

“ Tre-/3^eft “ 2 


(45) 


and the Gibbs state relative to the Hamiltonian Hes corrected by the Lamb-shift contribu¬ 
tion 

g-/3(Heff-hA2His) 

In general, £»st may differ significantly from the two Gibbs states above, their proximity or 
not being measured by the trace distance in terms of their Bloch expansion coefficients 

d{ga, gb) = ^Trv^ {ga - ghf = J (r? - r\f -k (rf - + (r^ - . 
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It can numerically be checked that the trace-distances between the Gibbs states (45) 
and (46) and the stationary state ^st are negligible under the same choice of parameters 
mentioned above in m- 

By the continuity properties of the relative entropy with respect to the trace-norm, it 
also follows that, under those physical conditions, the Definitionsandof internal entropy 
production practically coincide. In the following section we will first deal with the fate of 
the second law of thermodynamics in relation to the Redfield and the completely positive 
dynamics without assuming that the stationary state be of Gibbs form; that is we will 
use the expression (14) for the internal entropy production. Only when investigating the 
matter numerically we will fix the above mentioned physical conditions and consider the 
expression (0, taken advantage of the above mentioned proximity of the two definitions 
of entropy production. 


4 Complete positivity and entropy production 


Although, as discussed above, the stationary state of the micro-circuit reduced dynamics is 
explicitly time-dependent, gst{t) = R^{t) Qst R{t), in order to study the internal entropy pro¬ 
duction, we can consider the time-independent master equations and the stationary states 
^st- Indeed, the unitary connection between the physical and the rotated representation 
yields 


o-{g{t)) = -KTt\ht[g{t)](^logg{t)-loggst{t)'j 

= -KTt(h[g{t)](^\ogg{t)-loggsty] = a{g{t)) . 


(47) 


By means of the Bloch parametrization (37) of g{t) and of the matrix C in (38), one gets 

3 3 


i=l fi=0 


Oi 


ro = 1 


while the spectral representation of g{t) yields 
3 


^og g{t) = ^ I 1 +I log 


l+r{t)\ ^ 1 
2 




1 — r{t) 


where [r{t)]‘^ = Then, 


cr 


3 3 

{g{t)) = K 

i=l fi=0 


rAt) , 1 -|- r{t) rf , 1 -|- r®* 

~ 77 \ log ^ . , 

r{t) 1 — r{t) 1 —r®’'/ 


(48) 


where 
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Internal entropy prodnction: nnmerical analysis We first study the internal entropy 
production at t = 0 as a function of the initial state g. Three dimensional plots can 
be obtained by computing cr^g) in terms of two Bloch vector components: by setting rs = 0 
and plotting a{g) as a function of ri and r 2 , one sees that, in the case of the Redfield 
dynamics generated by there are regions where the entropy production is negative 

(see Fig. [^. None of these violations appear if the reduced dynamics is completely positive 
as that generated by C. 



Figure 2: (j{q) as a function of ri_ 2 , Redfield dynamics. 


Typically, the physical inconsistencies associated with the appearance of negative prob¬ 
abilities either in the spectrum of a time-evolving density matrix of the open quantum 
system (lack of positivity) or in the spectrum of an entangled state of the open quantum 
system and any finite level system (lack of complete positivity), manifest themselves at 
very short times. Quite different is the case with the second law of thermodynamics; its 
violations at time t = 0 are in fact not a negligible transient effect. Indeed, a numerical 
computation of (48) as a function of time, with physical parameters chosen as in |19) (see 


previous section), shows repeated violations of the second law also in the course of time. 
These violations occur for a same percentage of initial states as the violations at t = 0 and 
involve also states that initially have (t{p) > 0 and thus show no violations of the second 
law at t = 0 . 


As an example, take as initial state the one with Bloch components |g) = (1,0, 2r, r), 
r = —Ijy/b (cfr. ([^): it corresponds to an eigenstate of 0-3 and was studied in [19]. It 
exhibits an initial a{gt=o) > 0 followed by periodic violations of a{g{t)) > 0 (see Figure]^. 
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Figure 3: a{g{t)) as a function of time: initial state \g) = (1,0, 2r, r), r = —1/^5. 

In the case of a mixed state |^3) = (1,0, r 2 ,r 3 ), r 2 = 0.5 and rs = —0.4, that starts with 
(^{Qt=o) < 0, again periodic violations of a{g{t)) > 0 appear (see Figure Q. 


cr 

Ak 



Figure 4: a{g{t)) as a function of time: initial mixed state = (1, 0, 0.5,—0.4). 

In the graphs of both Figures and the blue line corresponds to the completely 
positive dynamics generated by £ that always yields a{g{t)) > 0 in agreement with the 
theory [8]. 

Remark 4. The observed violations of the second law of thermodynamics are not restricted 
to the specific choice of initial state and physical parameters considered in Figure]^ for sake 
of comparison with the experimental context devised in m- Indeed, the conflict between 
the non-complete positivity of the Redfield dynamics and the non-negative internal entropy 
production manifests itself across a whole range of parameters, namely for temperatnres 
between 0.0006 K and 0.06 K and ratios II/A between 0.1 to 10. In particular, viola¬ 
tions of the second law of thermodynamics at time t = 0 always occur, reaching 45% of 
the initial pure states at low temperature. Further, for every choice of temperature and 
pumping frequency, it is possible to find some initial state for which violations of the sec¬ 
ond law of thermodynamics occur repeatedly in time, becoming more and more typical for 
temperatures below 0.006 K. Violations of the second law of thermodynamics are therefore 
not exceptional, rather they are inherent to the non complete positivity of the considered 
Redfield dynamics. Whether such violations of the second law of thermodynamics are a 
feature of all non-completely positive dissipative dynamics is an open question; an answer 
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to it would demand either the proof that complete positivity is not only sufficient but also 
necessary to the non-negativity of the internal entropy production or devising an example 
of non-completely positive dissipative dynamics that does not conflict with thermodynamic 
expectations. Both tasks would require a stronger characterisation of the generators of pos¬ 
itive, but not completely positive dynamical maps, an issue which is still an open problem 
both mathematically and physically. 

4.1 Currents 

From an experimental point of view, due to the high time-resolution achieved by the present 
measnrement devices, discriminating the behavior of the internal entropy production in the 
Redfield and completely positive cases is in line of principle perfectly possible through a 
tomographic reconstruction of the time-evolving state. However, a more direct check of 
complete positivity in this model is possible by studying the asymptotic current supported 
by the micro-circuit. Indeed, the model was devised as to provide an experimental setting 
for probing the fate of electronic currents in open quantum micro-circuits. These specific 
characteristics can now be used to probe the complete positivity of the micro-circuit reduced 
dynamics; in fact, by numerical integration of the master equations made with the param¬ 
eters chosen in m, the temporal patterns of the cnrrents supported by the micro-circuit 
under the Redfield and completely positive dynamics appear different enough to allow for 
an experimental test. This fact may enable one to discriminate which one of the two pos¬ 
sible dynamics is actnally the best description for the true dynamics, thereby sorting out 
the fate of the second law of thermodynamics in the present model. 

The first of the two following figures shows the differences in the time-behaviour of the 
current supported by the two reduced dynamics starting with the initial state considered 
in Figure above. The second one refers to the fact that, despite the extremely high 
time-resolutions nowadays achievable experimentally, what one may hope to observe are 
not the pure oscillations in Figure but oscillations mediated over a few periods. Despite 
of this, as shown in Figure it looks possible to discriminate between the completely and 
non-completely positive patterns. 


I 

lo 



Figure 5: Time behaviour of the current supported by the micro-circuit; dashed lines represent the 
asymptotic valnes. 
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Figure 6: Time behaviour of the current supported by the micro-circuit mediated over two periods; 
dashed lines represent the asymptotic values. 


The possibility of discriminating between the two types of dynamics does not only hold 
for the specific values of the physical parameters considered above as chosen in uni; indeed, 
as shown in the following Figures and for temperatures T < 0.1 and rates Q/A in 
the range 0.3 — 10, the stationary states for the Redfield and completely positive dynamics 
support current behaviours different enough to be amenable to experimental tests. 



Figure 7: Asymptotic current as function of the temperature, with A = 0.005 and 12 = 2A and 
initial state as in Figure 
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Figure 8; Asymptotic current as function of il/A with A 
0.1 and initial state as in Figure|^ 


0.005, temperature such that /cT/A = 


5 Conclusions 

A typical argument against the request of complete positivity is that what really physically 
matters in the case of an open quantum system dynamics is the preservation of the posi¬ 
tivity of the eigenvalues of time-evolving density matrices, for they must be interpreted as 
probabilities. From this point of view, advocating the possible entanglement of the sys¬ 
tem of interest with an auxiliary inert system in order to justify the request of complete 
positivity looks like an abstract constraint. 

We have here showed that complete positivity cannot so easily be dismissed because 
of its connection with the second law of thermodynamics. Indeed, although it was known 
that complete positivity is sufficient to achieve non-negative internal entropy production, 
no concrete, experimentally appealing models have been proposed to check the connection 
between the structure of an open quantum dynamics and possible violations of the second 
law of thermodynamics. We have thus considered a model of open driven quantum micro- 
circuit which, in line of principle, can be experimentally studied, and showed that, if its 
dissipative dynamics is described by a Redheld dynamics, it would conflict with the principle 
of non-negative internal entropy production on a large variety of initial states, either at time 
f = 0 or repeatedly in the course of time. This latter fact is particularly significant as, in 
general, unphysical effects due to lack of complete positivity are thought to occur on very 
short time-scales (a similar situation in the case of non-Markovian effects is discussed in 

123 ]). 

The model was originally proposed to experimentally probe the electronic current sus¬ 
tained by the micro-circuit in weak interaction with its thermal environment: we have shown 
that looking at the temporal behaviour of the current may provide an experimental way to 
discriminate between the non-completely positive Redheld dynamics used in the model and 
the completely positive one that we have derived by adapting to a time-dependent Hamil¬ 
tonian the standard weak-coupling limits techniques. In turn, this would provide a check 
of the validity of the second law of thermodynamics at the quantum microscopic level. 

Indeed, the violations of the second law are not related to exceptional initial states or 
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values of the physical parameters of the model, rather they are typical and related to the 
structure of the Redfield dynamics which we have shown to be not even positivity-preserving. 
Whether the conflict with thermodynamics occur for general non-completely positive dy¬ 
namics is an open question: it is not known that complete positivity is also necessary to yield 
non-negative entropy production, nor there exist examples of only positivity-preserving dy¬ 
namics that do not violate the second law of thermodynamics. Indeed, on one hand, unlike 
for completely positive time-evolutions, there is no general characterisation of the genera¬ 
tors of only positive ones; on the other hand, a microscopic derivation of a non-completely 
positive dissipative dynamics is very likely to yield a not even positivity-preserving time- 
evolution. Hopefully, the findings of this manuscript will stimulate further research along 
these lines. 
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A Redfield type master equation 


In this appendix we discuss the derivation of the master equation (22): we follow the so- 
called projection technique 1211 [23 applied to the interaction representation. We thus set 


KsB{t) = , 


(49) 

(50) 


whence 

= XKtiRsBit)] , KtiRsBit)] = -i [KsBit ), R5B(t)] . (51) 

Let Qis be the bath Gibbs state at temperature T such that [Hq , g^] = 0; the following 
linear operators act as projectors on the states gsB{t) of the compound system S + B: 


P[R5B(t)] = (Trg ( Rsb 


QfS = R{t) (g) Qis 


= id — 


(52) 


with 


R{t) = TTB(RsBit)) = e^^^^^^TiBlgsBit)) e 




(53) 


giving the time-evolving density matrix of the open quantum system S in its own interaction 
representation. 


Let Kf ^ = 


o Ki o P, = P o Ki o 


rQP _ 


= Q o Kt o P and = Qo'Kf o > 


where o denotes the composition of maps. Then, (51) splits into the two coupled differential 
equations 


dP[R5g(t)] 

dt 

dQ[R5g(^)] 

dt 


XKfP o P[R5g(t)] + o Q[RsB{t)] 

XKf^ o P[R5B(i)] + AKf® o Q[RsBit)] ■ 


(54) 

(55) 
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The second equation is formally solved by 


Q[RsB{t)] = o + X ds Wgf o KfP o P[i 25 B(s)] 


(56) 


with the time-ordered solution to 


dW' 


QQ 

h* _ \ wQQ -mQQ 


, =XKY'^oWy , 

dt * ’ 


Ws.s = id 


With initial condition qsq = qb, from = 0 one gets 


(57) 


Q[RsB{t)] = A / ds Wgf o o P[R 5 b(s)] , 


(58) 


whence, inserting (58) into (54), 


'^El^k^M = XKfPoF[RsB{t)] + X‘^ f dsK^^ oWff oK‘^P oF[RsBis)] . (59) 

dt .In ’ 


The form of the interaction Hamiltonian in equation (21) and the fact that position operators 
have vanishing mean values with respect to Gibbs states yield TrB(^QBRsB{t)^ = 0; then, 
lP[RcSB(t)] = R{t) ® Qis implies 




KsB{t),Q o w9^ O Q Ksb(u) , R(u) (g) 


(60) 


The above equation depends on the history of the system state R(s) for all times 0 < s < t; 
in order to eliminate this dependence, one takes into account the weak-coupling hypothesis 
A <C 1 and looks at the dynamics as a function of a slow time parameter r = tX^. Firstly, 


by a change of integration variable s = t — u, (60) is recast as 

rt 


tQQ 


KsB{t),Q o „ O Q KsB{t - u) , Rt-u ® gp 


(61) 


Then, letting A —>■ 0, —)• id as the right hand side of (57) vanishes, and 


)oWi^foQ 


Ksb{s) , R{s) ® gp 


Ksb{s) , R{s) (g) gi3 


Ksb{s) , R{s) (g) gp 


The last equality follows from dhsi^gfiKsBis)^ = 0, as explained before. 

At this point, one usually sends the integration upper limit to +oo and, in R{t — u), 
replaces t — u = t/X^ — u with t; then, the equation (61) reads 


dR{t) 

dt 


r+oo 


-A' 


duTrg 


KsB{t), KsB{t - u), R{t) ® gfi 


(62) 
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By going back from the interaction picture to the Schrodinger one, the following master 
equation for g(t) is finally obtained, 


dg(0 

dt 


= ht[Q{t)] = -i +A2Nt[£.(t)] 


^+oo 

Moit)] = -J^ g{t)0Qp 


where 




(63) 

(64) 

(65) 


Using (21) and the thermal state 2-point functions (25) one finally obtains: 


«=1,3 


+ 


C*{u}n,u) g{t)e^^'^^[a^{t-u)], a^{t) | 


where 


C{ujn, u) = 2mw„Trf5 (gjj 


-iuHs piuHs 


( 66 ) 


(67) 


Since the couplings Xn do not depend on = 1,3 the explicit dependence on time t disap¬ 
pears. Indeed, let 

( cos nt 0 sin nt\ 

0 10 ( 68 ) 

— sinflf 0 cosQt/ 


be the matrix which implements the rotation (21): 

> ? = 1>3 . 

77=1,2,3 

Then, for generic 2x2 matrices A and B one finds 


(69) 




^=1,3 C=l,3 771,772=1,2,3 

= Y '^mri 2 i-'d^^mBar, 2 = Y ■ 


771.772=1,3 

Therefore, Nt[g(f)] becomes time-independent and equals 


J=l,3 


/— 1-00 

'2 / du- 


meit)] = -YY^ 

^=1,3 77 

+ C*icOn, u) [^ 7 (t) a^{-u) 


|c(wn, u) ct5(-m) e™^«ff^(t) 




}. 


(70) 


and, using (24) and (26), the master equation (63) reduces to (22). 
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Therefore, the generator Lt in ( |63[ ) becomes time-independent, too: Lt = L. In order to 
recast it in Lindblad form as in Q and Q, we first pass from the Pauli triple cri^ 2,3 to the 
rotated one, (Ti, 2 , 3 , in (36): 


j=i 


V = 


‘^eff 


0 A 17 
0 -17 A> 


(71) 


Then, ^ '5'3 and 




'^cos(a;efrt) — sin(a;efrt) 0^ 
^ Ul^{t)ak , ^eff(t) = I sin(a;efft) cos(a;efft) 0 

j,k=l 


(72) 


One can thus rewrite the term N[^>(t)] as follows: 


l^[^'(t)] = - X] [ duZjk{-u) ^G{u) dj,akQ{t) +G*{u) g{t) dk , aj , (73) 


j,k=l 


where, G{u) is as in (24). By taking into account that ^ = 1,3 the coefficients Zjk{t) can 
be regrouped into the following matrix by introducing the projection V = diag(l,0,1): 


/ cG+^sS ^cS-sG ^S\ 

* ^eff ^efF ^efF 


z{t) = vvn{t)v^Uesit) = 


QA 


sG- 

^cS 

^cC + 

0 

sS 


^^efF 




sG- 

^cS 

^cC + 

A 

sS 


<^eff 


‘^eB 


= cos 

(Ht), S 

= sin(Ht) 




OA 

A2 


G 


G 


(74) 


Finally, by separating the purely dissipative contribution D[^(t)] to N[£»(t)] from the one 
corresponding to a Lamb-shift Hamiltonian Hls, one gets the right hand side of (63) as 
follows: 


L[^(f)] 

Kjk 


— I 


H.s + X'Hls, git) +A 2 D[^)(f)] 


^ Kjk(dkgit)dj - , ^(t)|^ 

j,k=l 


^^'"du(G(T)Z,fc(-u) + G*iT)Zkji-u)) =K*kj , 


3 

where Hls = E Hjkdjdk, with 
j,k=l 


(75) 

(76) 

(77) 


Hjk — „. 


i^^“dr(G(r)Z,fc(-u) - G*iT)Zk,i-u))=H*k^ . 


(78) 


23 













In order to recast the action of L[^j(t)] as that of a 4 x 4 matrix —2£ on the Bloch vector 
(1, ri(t), r 2 (t), r 3 (t)) as in (37), one considers the linear action of L on the Pauli matrices 
<5 'i, 2,3 and on the identity do = 1: Ih[d^] = With ro(t) = 1 because of trace 

conservation, this gives 


L[e(t)l = I 


i=i 





(79) 

j=l 1^=0 




/ 0 0 0 




1^10 1^10 1^10 

1^10 


(80) 

^20 ^21 C.22 

C-23 



C = HeS +X^nLS = 

V-C30 -^31 ^32 C.33) 

The matrix C consists of an antisymmetric Hamiltonian contribution T-LeS + where 


'^eff — 


'^LS = 


/O 0 0 0\ 

0 0 0 
0 — Weff/2 0 0 


Vo 0 
/o 0 

0 0 

0 —'H12 

Vo —?7 i3 


-0 

0 

T~Li2 

0 

—'H 23 


0 / 

0 \ 
T~il3 
'H 23 
0 


plus a purely dissipative term 


V = 


/ 0 0 0 0 \ 
/Cio /Cii /C12 /Ci3 
/C20 ^12 ^22 ^23 
\JC30 JCi3 IC23 

ICn = K22 + K33 

/C22 = Kii + K33 

A ^33 = Kii + K22 


'Hi 2 = 2Xm{H2i) 
7ii3 = 2Zm{H3i) 
,7^23 = 2Xm(H32) 


'/Cio =Tm{K23) 
/C20 =Tm{K3i) 
/C30 =Zm{Ki2) 

'JC 12 =-TZe{Ki 2 ) 
/Ci3 =-7^e(^^l3) 
,Ai^23 = —T^e{K23) 


(81) 


(82) 


(83) 


(84) 


Using the expressions in (74), (77) and (78), the matrix entries explicitly read 

r*H-O0 


/Cio = 


A 

^eff Jo 


r*+oo 


dtt / dcu J(a;) sin(wn) ^ sin(a;effu) sin(Hu) + 

H- ( cos{uJesu) — 1 ) cos(Htt)') 

CUeff V / / 


10 


/C20 = - 


A r+°° 
^eff Jo 


r*+00 


/C30 = 


du 


du J dw J{uj) sm{oju) ^ — ^1 + cos(weff'^i)l^ sin(Hu)^ 

H-sin(a;eff^^) cos(Hu)') 

^efF ^ 

I dujJ(uj) sin(wu) (-^-sin(t<;effu) cos(Hu) + 

n V a; Vr 


+ 


+2-cos(a;eff^^) sin(Htt)') 

^efF ' 


(85) 


( 86 ) 


( 87 ) 


24 















/■+c« ^ 

dit / dcj J(w) cos(cjit) coth(—) (-sin(a;eff^i) sin(Ott) + 

Jo 2 VcJeflf 


/Cii = 2 


\ 

H—^ cos(ti;efr^i) cos(r2it) H-^ cos(Oit) j 


A2 


(88) 


/Cl2 = 


/Ci3 = 


/C 22 = 2 


/* + 00 /•+00 

-^ / du / da; J(a;) cos(a;tt) coth( —) sin(a;effM) cos(Oit) (89) 

■^eff -^0 io 2 

- / dtt / da;J(a;) cos(a;ii) coth(—) (-sin(a;efru) cos(r2u) + 

•^eff Jo Jo ^ ^ ^eff 

+ sm{Qu) ^ cos(a;efru) “ ) (90) 

r+00 r+00 

/ du j da; J(a;) cos(a;tt) coth(—) ( cos(a;efrtt) cos(i2ii) + 

Jo Jo 2 V 

n 


i I \ 

H-sin(a;efrti) sin(r2u) H- 7 ^ cos{^lu )) 


(91) 

(92) 


^ /*+oo r+oQ 1^^ . 

/C23 =-/ dw / do; J(a;) cos(a;u) coth( —) ( sin(a;effn) sin(llit) + 

i^efr Jo Jo 2 V 


0 JO 

H-^1 + cos(a;efru)^ cos{Qu)^ 

r*+oo r-\-oo 


(93) 


r+00 r+00 / 211 ^ + A ^ 

^33 ~ ^ / dn / do; J(a;) cos(a;tt) coth(—) (- 2 -cos(a;effu) cos(12it)+ 

Jo Jo 2 V a;^^ 


/o Jo 

n 


+2 -sin(a;effii) sin(12it)') ; 

a;cff- / 


^eff 

r + OO /* + CXD 


(94) 


JJl2 — 


+00 .+00 ^ ,/3a;, /2112 + A2 . ^ ^ x 

/ dii / do; j (a;) cos(a;u) cotn( — jl- 2 -sin(a;effw) cos(Sm)— 

JO Jo 2 V a;gg 


—2-cos(a;efrtt) sin(llii)') 

^efF ^ 


(95) 


r + OO /‘ + CXD 


'1^13 — 


- I du I do; J(a;) cos(a;u) coth(^) -sin(a;efr^i) cos(r2it) — (96) 

“^eff Jo Jo ^ ^^efF 

(97) 


— sin(lltt) + cos 


JJ 23 



sin(a;effii) sin(llr) + 


(98) 


In the main text, we have added a superscript “Red” to C,'Hls and T> in order to distinguish 
them from the analogous expressions pertaining to a completely positive dynamics which 
are obtained in the next Appendix. 
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B Completely positive master equation 


A physically consistent reduced dynamics can be obtained by a more careful treatment; it 
leads to a completely positive time-evolution, thus avoiding all the inconsistencies of the 
Redfield dynamics used in |T9j. We again use the projection technique but we follow the 
analysis of [8] without passing to the interaction representation. By repeating the arguments 


of the previous Appendix, one arrives at the following analog of equation (60): 
dg{t) 


dt 


— i [HeS ; i?(0] 

- / dnTi'B 
Jo 


Hssit), Q o o Q Hsb{u) , Q{t) (g) qb 


(99) 


As shown in [8] , a sounder strategy than the one that led to equation (62) in the previous 


Appendix consists firstly in formally integrating (99), yielding 
Q{t) = - 


-A^ du 
Jo Jo 


Trg 


Hs13{u),Q°'^uS , g{v) ® QB 


( 100 ) 


Secondly, in changing the double integral into 


du / 


'0 


f'dv f 

lo Jo 


Tre 

Trg 


Hsb {u),Q° ° Q Hsb {v), q{v)® qb 

Hsb{v + w) ,Qo o q Hsb{v) , q{v) ® QB 


and, finally, in going to the slow time-scale r = tX^, A << 1 where, using (65), one replaces 




Hsb{v) , q{v) ® QB 


by e 


I'ii’fHeff-l-BlB) 


Hsb{v) , Qiv) ( 8 ) qb 


so that the second integral with respect to drc becomes 
/• + 00 


N^[^>('f^)] = / dtue 
Jo 


This yields 


-■wHeft 


Tt, 


B 


Hsb{v + w ), Hsb{v) , Qv® Qb 


101 ) 


Q(t) = - A^ / duNt,[£i(u) 

Jo 


that solves the master equation: 

dg(t) 

dt 




H,fi,Q{t) +A2Nt[^>(t)] . 


( 102 ) 
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By proceeding as in the previous Appendix, one recasts (102) in the time-independent form 


dQ{t) 


dt 

n[git)] = -i 


= h[Q{t)] = -i Hes,Q(t) -hA^N[£»(t)] 


Hls , Q{t) + ^ Kjk(ak git) &j - -\^dj&k , 


(103) 

(104) 


j,k=l 


where Hls = k=i ^jkdjdk and, with respect to (77) and (78), the coefficients now read 

H,k = Yi ( 10 ^) 
K,k = du {g{u) Zkj{u) + G*{u) Zjkiu)) = . (106) 


/o 


Remark 5. In the Bloch representation, the generator L[^)(t)] corresponds to the action on 
the Bloch vector (1, ri(t), r 2 (t), r 3 (t)) of a 4 x 4 matrix C = Heg + YHls + YV, whose 
coefficients are all equal to the ones in (85)- (98) apart from JCiq, /C 12 , /C 13 and 7^23 which 
have opposite signs. 


The formal solutions to (103) are given by 


g{t) = + X^ dse®®® o N[£»s] 


(107) 


with e®®*^[£»] = exp(—itReff) Q exp(zt77eff)- 

On the slow time scale r = X^t, one rewrites 

g-iHeff Q e*^[£»] = g + j du o N o 1 






e'^^ig] 


(108) 


when A —>■ 0 the fast oscillations in the term in curly brackets average to zero. This allows 
one to replace that term by its ergodic average 


N = lim 4 ^ / ds ^ 

T— >--i-oo ‘2T y_7" 

which fulfils N o Heff = Heff o N, so that the resulting master equation is 

dg{t) 


dt 


= L[e(t)] = -i Reff + A^ Hls , Q{t) + A"D[^(t)] 


In the Bloch representation the action of N corresponds to that of the 4x4 matrix 

1 F — 

M= hm — dsUes{-s)AfUes{s) , 

T^+oo 21 Jq 


(109) 


( 110 ) 


( 111 ) 


with J\f the 4x4 matrix corresponding to N and Ues{s) the 4x4 matrix with 1, 0, 0, 0 in 
the first row and column and, in the rest, the 3x3 matrix in (72). Then, the action of L 
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can be represented by means of the 4x4 matrix C 
as in (81), while 


T-LeS + X^T-Lls + TX, where HeS is 


T-I-ls 



0 

0 

o\ 


( ° 

0 

0 

0 \ 

0 

0 

77i2 

0 

V = 

0 

^11 + K.22 

0 

0 

0 

—TLi2 

0 

0 


0 

0 

Alii -h JC 22 

0 

VO 

0 

0 

0^ 


V^30 

0 

0 

^33/ 


( 112 ) 
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